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FORWARD PION PRODUCTION IN p+p AND d+Au COLLISIONS AT STAR a 

G. RAKNESS, for the STAR Collaboration 
Penn State University /Brookhaven National Lab 



Measurements are reported of the production of high energy tt° mesons from the STAR ex- 
periment in p+p and d+Au collisions at ^/snn = 200 GeV and (77) = 4.00 (d beam direction). 
The inclusive yield agrees with perturbative QCD calculations in p+p collisions, but is re- 
duced in d+Au collisions. The azimuthal correlations of the forward 7r° with charged hadrons 
at midrapidity agree with PYTHIA in p+p collisions, but are suppressed in d+Au collisions. 
The results are consistent with the conjecture that the gluon density in nuclei is suppressed. 



1 Introduction 



In contrast to the nucleon, very little is known about the density of gluons in nucleiM J For protons, 
the gluon parton distribution function (PDF) is constrained primarily by scaling violations in 
deeply-inelastic lepton scattering (DIS) measured at the HERA collider. The data are accurately 
described by QCD evolution equations that allow the determination of the gluon PDF. It is 
found that the gluon PDF increases as the momentum fraction of the parton (xBj) decreases. In 
nuclei, the density of gluons per unit transverse area is expected to be larger than in nucleons. 
The interplay between large gluon densities and unitarity requirements on cross sections makes 
the nucleus a natural environment in which to establish if, and under which conditions, gluon 
saturation occurs. Quantifying gluon saturation is important because of the expectation that the 
nuclear matter created in central collisions of Au-Au nuclei evolves from an initial state produced 
by the collisions of the low- 2; Bj fields in each nucleus^ Fixed target nuclear DIS experiments 
have reported a suppression of the inclusive structure function normalized to proton DIS at 
\ow-xbj, but are limited in the kinematic range necessary to determine the nuclear gluon PDF. 

Many models exist which attempt to describe nuclear effects at low-x^j. Saturation models^ 
include a QCD based theory called the Color Glass Condensate (CGC), where gluon splitting is 
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balanced by gluon-gluon recombination to create dense macroscopic field J Another approach 
models quarks scattering coherently from multiple nucleons, leading to an effective shift in the 
XBj probedP Shadowing models suppress the nuclear gluon PDF in a standard factorization 
framework^ Parton recombination models modify the fragmentation of a quar k pa ssing through 
many gluonsP Other descriptions include factorization breaking in heavy nucleiP^With ^ ne data 
presently available, the mechanism by which low-x Bj gluon suppression occurs is unconstrained. 

Using factorization in a perturbative QCD (pQCD) framework, the PDF's and fragmen- 
tation functions (FF's) measured in electromagnetic processes can be used in the description 
of hadronic scattering processes. In p+p collisions at = 200 GeV, factorized leading twist 
pQCD calculations ha ve bee n shown to quantitatively describe inclusive ir° production over a 
broad rapidity wmdow ft^^ In pQCD, forward tt production in p+p collisions can be viewed as 
a probe of low-x^j gluons (g) in one proton with the valence quarks (q) of the other. Recently, 
the yield of forward negatively charged hadrons (h~) in the d-be am direction of d+Au collisions 
was found to be reduced when normalized to p+p collisions'*^ The reduction is especially sig- 
nificant since is ospi n effects are expected to reduce br production in p+p collisions, but not in 
d+Au collisions?^ The h~ suppression in d+Au has been interpreted as gluon suppression at 
low-x Bj ■ 

Insight on the particle production mechanism can be gained by analyzing the angular cor- 
relations of the forward hadron with a coincident hadron at midrapidityU"' A pQCD calcu- 
lation of the part on k inematics probed in inclusive forward hadroproduction reveals a broad 
distribution in xbJ^ Assuming collinear elastic parton (2 — > 2) scattering, the pseudorapidity 
(r/ = — ln[tan(#/2)]) of the second particle is correlated with XBj of the probed gluon, and a 
strong back-to-back azimuthal correlation is expected. In a saturation picture, the quark is 
expected to undergo multiple interactions through the dense gluon field, resulting in multiple 
recoil partons instead of a single recoil parto m' l -^ 



2 Experimental Results 

The Solenoidal Tracker At RHIC (STAR) is a multipurpose detector at BNL. One of its principal 

components is a time projection chamber used for tracking charged particles produced at < 

1.2. A forward tt° detector (FPD) comprising 7x7 matrices of 3.8 x 3.8 x 45 cm Pb-glass 

detectors was installed to detect high energy n° mesons with 3.3 < rj < 4.1. Data were collected 

over two years of RHIC operations at ^/snn = 200 GeV. In the 2002 run, p+p collisions were 

studied with a prototype FPD.-'" - In the 2003 run, p+p collisions were studied and exploratory 

measurements were performed with d+Au collision s. 

The cross sections for p+p — ► tt°+X at (77) = 3. 3^1 and (77) = 3. 8^*^ were reported previou sly. 

Preliminary results at (77) = 4.00 are compared with NLO pQCD calculations in Fig. ^ (left)M^ 

The largest ^-dependent systematic error at (77} = 4.00 is the energy calibration, resulting 

in 0.08 < (Jutt/ov < .12. The absolute angle uncertainty is the largest contribution to the 

normalization error?-^ The two curves use two sets of FF's, differing primarily in the gluon-to- 

pion FF, Dg. As pt decreases, the data make the transition from consistency with KKP to 
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consistency with Kretzer. This trend is also observed at midrapidityr^ At low px, the dominant 
contributions to 7r° production are qg and gg scatterings, making Dg the dominant FFp^l 

Nuclear effects on particle production are quantified by R^Am the ra tio of the inclusive 
yield of X in d+Au compared to p+p collisions normalized by the number of nucleon-nucleon 
collisions. This same measure was used at midrapidity, where it was found R£ Au > l^JJ The 
solid circles in Fig. ^ (right) are STAR preliminary data for R^ n at (77} = 4.00. On average, 
0.5 more photons per event are observed in d+Au collisions than in p+p collisions for events 
with > 30 GeV detected in the FPD. This leads to the largest pr-dependent systematic error 
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Figure 1: Inclusive tt° yield for p+p (left) and d+Au collisions normalized by p+p (right). The pion energy (E n ) 
is correlated with the transverse momentum (pr), as the FPD was at fixed values of pseudorapidity (77). The 
inner error bars are statistical, while the outer combine these with the E v - (pr-) dependent systematic errors, 
and are often smaller than the points. The curves (left) are NLO pQCD calculations evaluated at fixed r/, using 
different fragmentation functions. The x's and stars (right) are BRAHMS data for h~ production at smaller rj. 
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Figure 2: Coincidence probability as a function of azimuthal angle difference between the forward it meson and 
the leading charged particle at midrapidity. The left two distributions are simulations, using PYTHIA for p+p 
and HIJING for d+Au collisions. The right two plots are data. The curves are fits described in the text. 



which comes from the background correction to the d+Au yield. The ratio -R^Au a ^ = ^-00 is 
significantly smaller than i?i\ u at small er con sistent with the trend expected from models 
which suppress the nuclear gluon densit y! ' * * Linearly extrapolating R^j LU to 77 = 4, R^ u is 
systematically smaller, consistent with expectations of isospin suppression of p + p — » h' 

Exploratory measurements of the azimuthal correlations between a forward tt° and midra- 
pidity were completed for p+p and d+ Au collisions . The data analysis and detailed compar- 
isons to simula tion s were presented earlierP^l The left two plots in Fig. I^Jare simulations, using 
PYTHIA 6.22^21 f or p+p an d HIJING 1.381 24 for d+Au collisions. The right two plots in 
Fig. [2 are data. The leading charged particle (LCP) analysis selects the midrapidity track with 
the highest pr > 0.5GeV/c, and computes the azimuthal angle difference = cp^o — (pLCP 
for each event. The normalized distributions are fit with the sum of a constant and a Gaussian 
centered at ir. Correlations near A(f> = are not expected because of the large r\ separation 
between the tt and the LCP. The fit values are the areas under both the back-to-back peak (S) 
and the constant (B), and the width of the peak (175). 



PYTHIA reproduces most features of the p+p data. The back-to-back peak arises from 
2^2 scattering, resulting in forward and midrapidity jets that fragment into the 7r° and LCP, 
respectively. The width of the peak is smaller in PYTHIA than in the data, indicating that 
the momentum imbalance between the jets is too small in PYTHIA. This was also seen for 
back-to-back jets at the TevatronP^ 

The back-to-back peak is significantly smaller in d+Au collisions than in p+p, qualitatively 
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consistent with the coherent scattering-^ and CGO^ models. HIJING includes a model of 
shadowing for nuclear PDF's, and predicts a sizable back-to-back peak. This is not observed in 
the preliminary STAR data. Complete assessment of systematic errors is underway. 

In summary, inclusive forward ir° cross sections from p+p collisions at -y/s = 200 GeV are 
consistent with NLO pQCD calculations. Azimuthal correlations of the forward 7r° with midra- 
pidity are described by PYTHIA, which uses leading order pQCD with parton showers. The 
success of both of these calculations implies that forward 7r° production arises from partonic 
scattering in p+p collisions at yfs = 200 GeV. In d+Au collisions, the inclusive yield of forward 
7T° mesons is found to be significantly reduced compared to p+p. The r\ dependence of the 
reduction is consistent with models which suppress the gluon density in nuclei, in addition to 
exhibiting isospin effects at these kinematics. Exploratory studies suggest that the back-to- 
back peak of forward 7r° mesons with midrapidity is suppressed in d+Au relative to p+p 
collisions. This is qualitatively consistent with expectations that the particle production in a 
dense gluon medium differs from conventional leading-twist NLO pQCD expectations. More 
data for forward particle production and di-hadron correlations in d+Au collisions are required 
to reach a definitive conclusion about the existence of gluon saturation in the Au nucleus, or if 
an alternative description of the forward rapidity suppression is valid. A quantitative theoreti- 
cal understanding of the rapidity and px dependence of di-hadron correlations would facilitate 
experimental tests of a possible color glass condensate. 
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